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ABSTRACT 

We have updated our radially-resolved semi-analytic models of galaxy formation to 
track the radial distribution of stars, metals, atomic and molecular gas in galactic 
disks. The models are run on dark matter halo merger trees from the Millennium and 
Millennium II simulations using the physical recipes outlined in Fu ct al. (2010) and 
Guo et al. (2011), with 3 main changes: (1) We adopt a simple star formation law 
where the star formation rate surface density is proportional to the surface density of 
molecular gas (Ssfr. °c Sjj 2 ). (2) We inject a fraction of the heavy elements produced 
by supcrnovae directly into the halo, instead of first mixing them with the cold gas 
in the disk. (3) We include radial gas inflows in disks using a model of the form 
^inflow = or. 

The average surface density profiles of molecular gas in L* galaxies place strong 
constraints on inflow velocities, favouring models where inflow ~ 7 km/s at a galac- 
tocentric radius of 10 kpc. The radial inflow model has little influence on the gas 
and stellar metallicity gradients in the outer regions of galaxies. Gas-phase metallicity 
gradients are affected much more strongly by the fraction of metals that are directly 
injected into the halo gas, rather than mixed with the interstellar cold gas of the 
galaxy. Metals ejected out of the galaxy at early epochs result in late infall of pre- 
enriched gas and flatter present-day gas-phase metallicity gradients. A prescription in 
which 80 % of all the metals produced by stars are injected into the halo gas provides 
the best fit to the relatively flat observed metallicity gradients of galaxies with stellar 
masses greater than 10 10 M©. Such a prescription also results in a good fit to the rela- 
tion between gas-phase metallicity and specific star formation rate in the outer parts 
of galactic disks. We examine the correlation between gas-phase metallicity gradient 
and a number of global galaxy properties, finding that it is most strongly correlated 
with the bulge-to-disk ratio of the galaxy. This is because gas is consumed when the 
bulge forms during the galaxy merger, and the gas-phase metallicity gradient is then 
set by newly-accreted gas. These model predictions appear to be in good agreement 
with observations from Moran et al. (2012). 

Key words: galaxies: evolution - galaxies: formation - stars: formation - galaxies: 
ISM - ISM: atoms - ISM: molecules 



1 INTRODUCTION 

In recent years, observations of radially resolved profiles of 
atomic gas, molecular gas and star formation in represen- 
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tative samples nearby galaxies (e.g THINGS for HI pro- 
files, Walter et al. 2008; HERACLES for H 2 profiles, Leroy 
et al. 2009) have motivated galaxy formation theorists to 
include the physics of the atomic-to-molecular gas transi- 
tion in their models. Some of these models make predic- 
tions for the global atomic and molecular gas content of 
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galaxies (e.g. Obreschkow et al. 2009; Lagos et al. 2011), 
while others attempt to model the radial structure of the 
gas in more detail (e.g. Robertson & Kravtsov 2008; Fu et 
al. 2010; Power, Baugh & Lacey 2010; Feldmann, Hernandez 
& Gnedin 2012). Radial abundance gradients also place im- 
portant constraints on disk formation models, particularly 
on how supernova feedback (SN feedback) processes eject 
metals into the surrounding gas and how these metals are 
mixed throughout the disk as the galaxy evolves. A success- 
ful model of disk galaxy formation should be able to repro- 
duce the metallicity gradients in galaxies along with radial 
density profiles of old stars, young stars and gas. 

Observational data on metal abundance gradients in 
galaxies remains rather sparse. Many studies focus only on 
one galaxy or on a handful of galaxies (e.g Rudolph et 
al. 2006 for the Milky Way; Magrini et al. 2007 for M33; 
Bresolin et al. 2009 for M83). For many years, the Zaritsky, 
Kennicutt & Huchra (1994) study of gas-phase metallicity 
gradients derived for 159 HII regions in 14 spiral galaxies 
in combination with published data for another 25 galax- 
ies, has remained the standard reference in the field. More 
recently, Moustakas et al. (2010) published metallicity gra- 
dients for 65 galaxies from the SINGS survey and Moran 
et al. (2012; hereafter Moranl2) analyzed gradients for 174 
galaxies with atomic and molecular gas mass measurements 
from the GASS and COLD GASS surveys (Catinella et al. 
2010; Saintonge et al. 2011a). These new observations make 
it possible to carry out a statistical comparison with model 
predictions for the first time. 

There is a long history of models of the radial structure 
and properties of disk galaxies in the literature, beginning 
with Tinsley & Larson (1978) who implemented a model 
for chemical evolution into the dynamical collapse calcula- 
tions for gas clouds with rotation and axial symmetry in- 
troduced by Larson (1976). Similar modelling efforts were 
later undertaken by many others (e.g. Matteucci & Fran- 
cois 1989; Kauffmann 1996; Chiappini, Matteucci & Grat- 
ton 1997; Dalcanton, Spergel & Summers 1997; Avila- Reese, 
Firmani, & Hernandez 1998; van den Bosch 1998; Prantzos 
1999; Dutton et al. 2007; Fu et al. 2009; Yin et al. 2009; 
Cook et al. 2010; Fu et al. 2010). The models can be ar- 
ranged in increasing order of complexity, from those that 
use simple parameterized formulae to describe the infall of 
gas onto the disk, to those that are embedded within high 
resolution N-body simulations of structure formation in a 
ACDM cosmology and where gas infall rates are determined 
by explicitly following the growth of each dark matter halo 
in the simulation and then calculating the rate at which gas 
will cool as a function of time. 

In order to model the radial distributions of stars and 
gas in the disk, processes such as star formation, SN feedback 
and chemical enrichment need to be treated in a radially re- 
solved fashion. In previous work (Fu et al. 2010, hereafter 
FulO), we developed models in which each galaxy disk was 
divided into a series of concentric rings. We adopted two 
different prescriptions to partition the cold gas into atomic 
and molecular components: one is adapted from the models 
of Krumholz, McKee & Tumlinson (2009), in which the H 2 
fraction is parameterized as a function of local cold gas sur- 
face density and gas-phase metallicity. The other prescrip- 
tion is empirically-based and proposes that the H2 fraction is 
related to the interstellar pressure (Elmegreen 1989 & 1993; 



Blitz & Rosolowsky 2004 & 2006). The FulO semi-analytic 
model is based on the version of the L-Galaxies code de- 
scribed in De Lucia & Blaizot (2007) (hereafter DLB07) 
implemented on the halo merger trees of the Millennium 
Simulation outputs (Springel et al. 2005). As discussed in 
the paper, the model reproduces the radial profiles of stellar 
mass, HI, H2, and star formation in L, galaxies. However, 
the model produces a stellar mass function that is too steep 
at the faint end and the gas fractions of low mass galaxies 
do not agree with observations. 

The analysis presented in this paper is based on the new 
version of the L-Galaxies code described in Guo et al. (2011) 
(hereafter Guoll), which runs on the halo merger trees of 
both the Millennium Simulation (Springel et al. 2005) and 
on the 125 times higher resolution Millennium II Simula- 
tion (Boylan-Kolchin et al. 2009). This allows us to study 
the formation and evolution of galaxies with masses ranging 
from those of dwarfs to the most massive cD galaxies. 

Our paper is organized as follows. In Section 2, we 
briefly describe the N-body simulation, and the changes to 
the semi-analytic model with respect to FulO and Guoll. 
We discuss how we normalize the free parameters in our 
models using gas profile data from Leroy et al. (2008) for 
a small sample of nearby disks, along with the stellar and 
gas mass functions of nearby galaxies at z = 0. In Section 
3, we present stellar, HI and H2 mass functions at z = 
and compare them to observations. In Section 4, we ana- 
lyze how the star formation rate surface density profiles of 
galaxies depend on stellar mass. In Section 5, we study how 
the radial profiles of gas-phase metallicity depend on stellar 
mass, bulge-to-disk ratio, stellar surface density and specific 
star formation rate. We also look at the relations between 
metallicity and star formation in the inner and outer disks of 
galaxies. The model results presented in Sections 4 and 5 are 
compared to observational results from Moranl2. Finally, in 
Section 6, we summarize our results and discuss avenues for 
future work. 



2 THE MODELS 

2.1 The N-body simulations 

The galaxy formation models are run on two very large 
N-Body simulations: the Millennium Simulation (hereafter 
MS, Springel et al. 2005) and Millennium-II Simulation 
(hereafter MS-II, Boylan-Kolchin et al. 2009). Both simu- 
lations adopt a ACDM cosmogony with parameters Q,a = 
0.75, Q m = 0.25, fi baryon = 0.045, cr 8 = 0.9 and h = 0.73, 
and they track 2160 3 « 10 10 particles from z = 127 to z = 0. 
The main difference between the two simulations is the res- 
olution. The particle mass in MS is 8.6 x 10 8 M Q h' 1 and 
the periodic box is 500 Mpc h~ on a side; the particle mass 
in MS-II is 6.8 x 1O 6 M h' 1 and the size of the periodic box 
is 100 Mpc h~ , which means the resolution of MS-II is 125 
times higher than MS. Since the smallest halo traced in both 
simulations has 20 particles, MS can be used to study the 
formation and evolution for Milky Way-sized or larger galax- 
ies, while MS-II is more appropriate to study dwarf galax- 
ies. For our models of the cold gas components in galaxies, 
the MS is suitable to study atomic and molecular gas for 
disk galaxies with stellar mass M* > 1O 1O M and MS-II 
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is useful to study the gas components in present-day dwarf 
galaxies and in the high-redshift progenitors of present-day 
L* galaxies. Another difference between the two simulations 
is the number of output snapshots: the last 60 snapshots 
for MS and MS-II are identical, but MS-II has more snap- 
shots at very high redshift to increase the time resolution in 
early- forming first structures. 



2.3.1 The radial distribution of gas cooling onto the 
galaxy disks 

We adopt the same prescription described in FulO, in which 
newly infalling cold gas at each time step has exponential 
form with scale length nnfaii = J"d = (A/v2) r v - a (Mo, Mao 
& White 1998). Gas is simply added to each radial ring such 
that it follows this profile. 



2.2 The galaxy formation models 

The reader is referred to Section 3 of Guoll for a detailed 
discussion of how the various physical processes in the model 
are calculated, including re-ionization, gas infall and cool- 
ing, star formation and metal production, SN feedback, ram- 
pressure stripping of gas in satellite galaxies, tidal disrup- 
tion of satellites, galaxy mergers, bulge formation, black hole 
growth and AGN feedback. 

The cosmological parameters for both MS and MS-II 
and the semi-analytic models in DLB07 and Guoll are from 
WMAP1 (Wilkinson Microwave Anisotropy Probe) first- 
year results (Spergel et al. 2003). To update the cosmology, 
Guo et al. (2013, hereafter Guol3) adopt the technique de- 
scribed in Angulo & White (2010) to rescale the growth 
of structure to be appropriate for WMAP7 parameters 
(fi A = 0.728, n m = 0.272, f2 baryon = 0.0454, ct 8 = 0.807 
and h = 0.704; Komatsu et al. 2011). We make use of these 
rescaled models in this paper, but we note that differences 
for the quantities explored in this paper are negligible com- 
pared to systematic uncertainties in our treatment of the 
various gas-physical processes. 

The main changes made in this paper with respect to 
Guoll are the following: (1) We include the methodology to 
track gas and stellar components in radial rings described in 
FulO. (2) We include the prescriptions for the transition be- 
tween atomic gas and molecular gas described in FulO. (3) 
We abandon the two-regime star formation model in FulO 
and adopt a simpler prescription where the star formation 
surface density is always proportional to the molecular gas 
surface density (e.g. Leroy et al. 2008; Bigiel et al. 2008; 
Schruba et al. 2011). (4) We allow a fraction of heavy ele- 
ments to be mixed directly with the hot gas in the halo. We 
describe these modifications in more detail below. 



2.3 Treatment of the radial distribution of gas in 
disks 

In this paper, we divide each disk into 12 rings instead of 
30 rings as in FulO. This is done to decrease the computa- 
tional memory requirements. The ring radii are given by the 
geometric series 



r% = 0.44 x 1.5* [h~\pc] (i = 1, 2. ..12) 



(1) 



The radius for the innermost ring in Eq. {JXJ) is about 0.9 kpc, 
and the radius for the outermost ring is about 80 kpc. As we 
mentioned in FulO, the radial profiles are insensitive to the 
precise scheme, if the adopted number of rings is sufficiently 
large. This sub-division of the disk in Eq. {TJ is sufficient to 
study both small and large galaxies. 



2.3.2 Mergers, bulge formation and the radial distribution 
of bulge stars 

The DLB07 model does not include a prescription for bulge 
sizes and the FulO model did not attempt to calculate the 
radial distribution of bulge stars. The Guoll models include 
a method for calculating the sizes of bulges produced during 
galaxy mergers and by disc instabilities. 

The size of the bulge after a merger is calculated using 
the following equation 



C- 



G h C h a- 



''i 



n + r 2 



(2) 



where C parametrizes the binding energy of the galaxy and 
a parametrizes the effective interaction energy deposited in 
the stellar components. In Guoll, C = 0.5 and a — 0.5 are 
adopted for galaxy mergers. m now and r ncw in Eq. ([2]) are 
the mass and half mass radius of the newly formed bulge. 
For major mergers, in which the ratio of the total (stel- 
lar +gas) mass of the satellite to the central galaxy is larger 
than 0.3, both disks are destroyed. The stars in both pro- 
genitor galaxies and the stars formed in the merger-induced 
starburst become bulge stars. In this case, mi is sum of 
the stellar mass of the two progenitor galaxies and m2 is the 
sum of the stellar mass converted from cold gas in the major 
merger starburst (see later) ; r\ and r-z are the corresponding 
half-mass radii. 

For a minor merger (m sa t/m cen < 0.3), the stellar com- 
ponent of the smaller satellite galaxy is added to the bulge 
of the central galaxy, while the cold gas component of the 
satellite is added to the disk of the central galaxy, mi in 
Eq. ([2| is the bulge mass of the central galaxy before the 
merger and is the stellar mass (bulge+disk) of the satel- 
lite galaxy; r\ and r2 are the corresponding half mass radii. 

The main new prescription implemented in this paper 
concerns the radial profile of gas added to the central galaxy 
from the satellite before the starburst occurs. The FulO code 
directly superposes the gas radial profiles of two gas disks 
together and then processes the starburst. In this paper, we 
treat the accreted gas from the satellite in the same way as 
gas cooling from the halo; the scale length of the accreted 
gas is determined from the spin parameter and virial radius 
of the halo of the central galaxy using the equation ri n f a ii = 
(A/V2>vir. 

The DLB07 and Guoll codes use the "collisional star- 
burst model" introduced by Somerville, Primack & Faber 
(2001). In both major and minor mergers, the mass of stars 
formed in the starburst is m* = eb U rstm gas , with eb ur st given 
by the equation 



eburst = Eburst (m sa t /m C en) C 

with Qburst = 0.7 and /3burst = 0.56. 



(3) 
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In our model, we simply calculate the surface mass den- 
sity of stars formed in the burst in ring number i as 



(4) 



using Eq. Q to define eburst- F° r starbursts occurring dur- 
ing major mergers, the stars formed in all the rings become 
bulge stars. 

This is admittedly not a very realistic description of 
what happens in an actual merger-induced starburst, but 
since we mainly study quiescent disk galaxies in this paper, 
this will not concern us. The other update we make with 
respect to the Guoll models during mergers is to calculate 
the half mass radii directly from the stellar and gas profiles, 
rather than to assume an exponential form. 

Bulges also form through secular evolution. When 
«max < \J Grn» id i sk /r» id isk, the disk is considered unstable. 
A mass <5m* from the inner part of the disk is transferred to 
the bulge 



<5m* = m* 



r*,disk^n 

G 



(5) 



where ra^disk and r*,disk are the mass and scale length of the 
stellar disk. If this material forms a new bulge, it is assumed 
to have a half-mass radius Sr equal to the outer radius of 
the region. If there is already a bulge in the galaxy, then the 
resulting bulge is calculated using 



Or ^ucw 
.5 



0.5— +0.5— +2.0 — (6) 

' new 

where mi and n are the mass and half mass radius of the 
existing bulge; m,2 = 5m, and ri = Sr are the mass and 
radius of the transferred region. Once again we have up- 
dated this prescription by calculating the radius Sr directly 
from the stellar radial distribution of the disk rather than 
by assuming an exponential stellar disk profile as in Guoll. 

As in Guoll, we assume a Jaffe profile (Jaffe 1983) for 
the bulge of the form 



Pbuigo (r) 



m*bulgc 



n, 



r 



(7) 



in which mb u i ge is the bulge mass and r^ is the half mass 
radius of the galaxy bulge. The stellar mass in a given radial 
ring is the combined mass of disk and bulge stars in that 
ring: 



m» = m„,disk + 47r / pbuige (r)r dr 



■ m, disk + mbule 



1 + 



- 1 + 



Tout 
^b 



(8) 



where m*,disk is the stellar mass of the disk component in 
the ring, and r- m , r out are the inner and outer radii of the 
radial ring. We note that the stellar surface density profiles 
presented in the next sections include both bulge and disk 
components. 



2.3.3 Radial gas inflows in disks 

In FulO, we showed that a star formation law of form 
Ssfr oc Eh 2 was problematic in that it led to cold gas sur- 
face density profiles that were too shallow to match obser- 
vations (see Sec. 3.4 & Fig. 3 in FulO). In that paper, we 



introduced the following two tweaks to solve the problem: 
(i) we adopted a star formation law in regions of the disk 
with /h 2 < 0.5 of the form Esfr oc £ gaa . (ii) we assumed 
that the SN reheating efficiency was inversely proportional 
to gas surface density (Am re heat oc Am*/E gas ). The second 
assumption increases the gas consumption timescale in the 
inner disk. 

In this paper, we will introduce a radial gas inflow pre- 
scription as a more realistic solution, i.e the gas from the 
outer disk flows inwards towards the inner disk. As we will 
now show, a plausible inflow prescription combined with a 
simpler star formation law in which the star formation rate 
surface density is always proportional to the molecular gas 
surface density, produces gas surface density profiles that 
match observations without need for a radially-dependent 
SN feedback efficiency. 

There are a number of physical mechanisms that can 
drive radial inflows of gas in the disk. The ones that are in- 
voked most frequently involve gravitational interaction be- 
tween gas in the disk and non-axisymmetric stellar struc- 
tures such as bars and spiral structures (Kalnajs 1972). Sim- 
ple physical considerations yield estimates of flow velocities 
ranging from 0.1 to a few km/s (Lacey & Fall 1985; Bertin & 
Lin 1996) . Attempts to measure radial flow rates in galaxies 
have been confined to very small samples and have yielded 
inconclusive results (e.g. Wong, Blitz & Bosma 2004; Haan 
et al. 2009; Zhang & Buta 2012). Difficulties arise from the 
fact that flow patterns in individual galaxies are often irreg- 
ular and disks are frequently not axisymmetric. In addition, 
different methods for estimating mass inflow rates yield dis- 
crepant results. 

As a result, modellers generally resort to simple param- 
eterized inflow prescriptions (e.g Lacey & Fall 1985; Porti- 
nari & Chiosi 2000; Schonrich & Binney 2009; Spitoni & 
Matteucci 2011). In this work, we have tested a number of 
inflow prescriptions. We first tried the simplest prescription 
in which the radial gas inflow velocity is constant for the 
whole disk (model (a) in Lacey & Fall 1985), but we found 
that the inflow velocities required to move enough gas from 
the outer disk to the inner disk to compensate the depletion 
by star formation, lead to an unacceptably large pile-up of 
gas in the very inner disk. A prescription in which the rate 
of change of the angular momentum is proportional to the 
angular momentum yields results that agree best with ob- 
servational data, i.e 



dL„ 



dt 



(0) 



where C is a constant for all galaxies. Because L gas = 
m gas r gas u c ir, we have that 



a v r 



(10) 



The radial inflow velocity of the gas is thus proportional 
to the galactocentric radius r. This is equivalent to model 
(b) in Lacey & Fall (1985). Spitoni & Matteucci (2011) also 
invoke a radial flow prescription of this form to reproduce 
the observed stellar metallicity gradient in the Milky Way. 

The constant a v in Eq. (|10[) is an adjustable model pa- 
rameter. As discussed later, a v — 0.7 km s _1 kpc -1 is re- 
quired to fit observations. This results in negligible flow ve- 
locities in the inner disk and flow speeds of several km s _1 
in the outer disk (see Fig. [1} . 
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Figure 1. The relation between galactocentric radius r and radial 
gas inflow velocity Uj n fl ow from Eq. IIIOH (a v = 0.7 km s~ x kpc -1 
is adopted). 



for the radial inflow of gas has a much weaker influence on 
metallicity gradients than the fraction of the metals pro- 
duced by SNe we choose to put into the hot gas rather than 
the cold interstellar medium. We see from Fig. [3] that inclu- 
sion of radial gas inflow yields stellar metallicity gradients 
that are slightly steeper in the region of the disk interior to 
3 kpc, but otherwise the predicted metallicity gradients are 
virtually indistinguishable in the two cases. We note that 
the Spitoni & Matteucci (2011) claim that radial inflows 
are required to explain the observed metallicity gradients in 
the Milky Way, but in their models accreted gas always has 
primordial metallicity. In our model, infalling gas from the 
halo has already been pre-enriched by the ejection of metals 
during the early formation phase of the galaxy. 



In the models, we calculate the gas inflow velocity ac- 
cording to Eq. (|10|l after processing gas cooling and infall in 
each time step and derive the new radius r — r — Wi n fl ow Ai 
for the gas in each ring, where At is the length of the time 
step. Proceeding from the inner ring to the outer ring, we 
move the gas to its new radius; metals associated with the 
cold gas move inward together with the gas. 

In Fig. [2j we show the best-fit radial surface density 
profiles of stars, atomic gas, molecular gas, and total cold 
(atomic+molecular) gas for disk galaxies in haloes with 
circular velocities comparable to that of the Milky Way 
(200 < v C i rc < 235 km/s) at z = 0. The four panels on 
the left show the results if radial gas inflows are not in- 
cluded, and the four right panels show the results of our 
inflow model. Without inflows, both the molecular gas and 
average total cold gas radial surface density profiles are too 
fiat in comparison to the data of Leroy et al. (2008). The 
right panels show that the inclusion of radial flows leads to 
gas profiles that are in much better agreement with observa- 
tions. We note that the stellar mass profiles are insensitive 
to the radial flow prescription; this is because the majority 
of the mass in the inner galaxy forms at high redshifts from 
gas cooling in smaller, denser haloes. The gas radial profiles 
are set by gas that cools at late times. Because the timescale 
over which H2 is transformed into stars is ~ 2 Gyr (see next 
section), the H2 surface densities in the inner disk quickly 
drop to low values as a result of star formation and SN feed- 
back, unless gas is transported from the outer disk to the 
inner regions of the galaxy. 

We note that we have tuned the parameters of the ra- 
dial inflow model (in particular the constant a v in Eq. I10[) 
so that the gas radial profiles in Milky Way-type disks are 
reproduced. One might ask whether metallicity profiles pro- 
vide an additional check on the inflow model. In Fig. |3j 
we show the mean radial gas-phase metallicity profile (left 
panel) and stellar metallicity profile (right panel) for the 
same set of Milky Way type galaxies at 2 = 0. Green curves 
show results with radial inflow and blue curves show models 
without radial gas inflow. In the left panel, the black dots 
represent oxygen metallicity measurements from HII regions 
in the Milky Way by Rudolph et al. (2006) and black dashed 
line is the corresponding linear fit to the data between and 
15 kpc. In the right panel, the black dots represent the stel- 
lar oxygen measurements gradients derived from Cepheids 
by Andrievsky et al. (2002a,b,c;2004) and Luck et al. (2003). 
Red dots with error bars are the corresponding mean values. 

As we will discuss in the next section, the prescription 



2.4 Star formation prescription 

In this paper, we adopt a simple law in which the star forma- 
tion surface density is proportional to the H2 surface density: 
Ssfr oc Eh 2 (e.g Leroy et al. 2008; Bigiel, Leroy & Walter 
2011; Schruba et al. 2011; Leroy et al. 2013). We adopt a con- 
stant molecular gas consumption timescale for all galaxies at 
all redshifts. There is evidence that molecular gas consump- 
tion times may be longer in massive, early-type galaxies than 
in normal, present-day star-forming spirals (Saintonge et al. 
2011b) and shorter in starburst galaxies (Genzel et al. 2010), 
but we will neglect these complications, because it is still un- 
clear what the underlying physical cause of these effects are. 

The star formation law adopted in the model can be 
simply written as 

Ssfr = Qh 2 Sh 2 (11) 

in which the H2 star formation efficiency an 2 = 5.3 x 
10~ 10 yr -1 is slightly adapted to yield a better fit to the 
results stellar mass function at z = 0. We note that the 
molecular gas surface density Eh 2 in Eq. (|lll) does not in- 
clude the helium component. 

2.5 Prescriptions for the transition between 
atomic and molecular gas in the ISM 

In FulO, we implemented two prescriptions for the transi- 
tion between atomic and molecular gas in the interstellar 
medium. One is from Krumholz et al. (2009) (H2 fraction 
prescription 1), who calculate an equilibrium H2 fraction 
for a spherical cloud with given dust content surrounded by 
a photo-dissociating UV field. In this prescription, the H2 
fraction is primarily a function of local cold gas surface den- 
sity and metallicity. The second prescription originates from 
Elmegreen (1989, 1993) and Blitz & Rosolowsky (2006), in 
which the H2 fraction is a function of the pressure in the 
ISM. We adopt the approximation by Obreschkow & Rawl- 
ings (2009), in which the local ISM pressure is expressed as 
a function of gas and stellar surface density in the radial 
ring. More details can be found in Section 3.2 of FulO. 

In this paper, we make a number of adaptations to the 
Krumholz et al. prescription. We update the prescription 
using the recent fitting equations in McKee & Krumholz 
(2010) with the molecular gas fraction /h 2 given by 
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Figure 2. The radial surface density profiles of stars, atomic gas, molecular gas, and total cold gas for disk galaxies similar to the Milky 
Way. Results from the model with H2 fraction prescription 1 are shown at at z = 0. The two panels show the model results with (right 
panel) and without (left panel) radial gas inflow. The light blue curves with error bars are taken from the data of Leroy et al. (2008) for 
disk galaxies in the range of circular velocities 200kms- 1 < v ciy . < 235kms- 1 (NGC 0628, NGC 3184, NGC 5194, NGC 3521). The red 
solid and dotted curves are the mean and the median values from model results, and the black dashed curves are the ±1<t deviations 
about the mean values. 
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Figure 3. The mean radial mctallicity profile of gas (left panel) and stars (right panel) for a Milky Way-sized galaxy at z = from 
models with (green curve) and without (blue curve) radial gas inflow. In the left panel, the black dots are the oxygen metallicity from 
HII regions in the Milky Way (Rudolph et al. 2006) and black dashed line is the corresponding linear fit to the data between and 15 
kpc. In the right panel, the black dots represent the stellar oxygen gradients from Cepheids (Andrievsky et al. 2002a, b,c;2004; Luck et 
al. 2003) and the red dots with error bars are corresponding mean values. 



for s < 2 and /h 2 = for s > 2. The s in Eq. Q is defined 

as 

ln(l + 0.6 X + 0.01 X 2 ) , 1Q , 
0.6t c 

in which X = 3.1 (l + 3.1Z' 365 ) /4.1 and r c = 
0.66 (£ comp /M pc" 2 ) Z'. Note that Z' = Z gSLS /Z Q is the 
gas-phase metallicity relative to the solar value, and E comp 
is the gas surface density of the giant gas cloud. Since the 
gas surface density in the model is the azimuthally averaged 
value in each concentric ring, a clumping factor Cf is intro- 
duced to consider the difference between E CO mp and E gas 

Ecomp = CfSg as (14) 

One problem with the Krumholz et al. prescription is 
that it can easily yield non-convergent results. The reason 
for this is that the value of /h 2 is very sensitive to the ex- 



act value of the gas-phase metallicity when the metallic- 
ity is well below the solar value (see Figure 5 in McKee & 
Krumholz 2010). Molecular clouds do not form until metals 
have been produced, and metal production is dependent on 
the formation of stars, which can only take place in molec- 
ular clouds, so the conditions for star formation to occur in 
the first generation of haloes to form in the simulation are 
subject to considerable uncertainty. Moreover, at the resolu- 
tion limit of the simulation, the assembly histories of haloes, 
and hence the star formation histories of galaxies embedded 
within them, cannot be tracked accurately. It thus becomes 
very difficult to get results that are consistent for galaxies 
forming in haloes of the same mass in simulations of differ- 
ent resolutions (e.g. MS and MS-II). Similar problems were 
reported by Kuhlen et al. (2012) in their attempts to model 
dwarf galaxies with Krumholz et al. prescription. 

In the top panels of Fig. [4l we plot the relation between 
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Figure 4. The relation between gas phase metallicity and galaxy halo mass at redshifts 3, 2, 1, and 0. The red curves and green curves 
are the mean values derived for the MS and the MS-II simulations, respectively. The four top panels are from the Krumholz et al. H2 
prescription with fixed clumping factor Cf = 1.5, and the four bottom panels are from the Krumholz et al. H2 prescription with variable 
clumping factor dependent on gas-phase metallicity (see text). 



the gas-phase metallicity of a central galaxy and the mass 
of its parent halo at different redshifts. The red curves and 
green curves show results for the MS and MS-fl simulations, 
respectively. We have adopted Cf = 1.5 for all galaxies at all 
redshifts, as in FulO. As can be seen, at fixed halo mass, the 
gas-phase metallicity is higher in the MS-II run than in the 
MS run, particularly for low mass haloes at high redshifts. 
Since the halo resolution of MS-II is 125 times higher than 
MS, a larger fraction of the formation and chemical enrich- 
ment history of low mass haloes lies below the resolution 
limit for MS than MS-II. Because the Krumholz et al. pre- 
scription is so sensitive at low metallicities, the resolution 
discrepancies are enhanced for high redshift galaxies in low 
mass haloes. 

The question arises as to why stars are able to form at 
all in low-metallicity galaxies. For example, 1ZW18 is the 
most metal-poor galaxy known, yet it is a starburst system 
(e.g. Martin 1996). The main reason is that the gas in low- 
metallicity dwarf galaxies is clumpy on scales of ~100 pc 
(Lo, Sargent & Young 1993; Stil & Israel 2002). To account 
for this, Kuhlen et al. (2012) adopt a clumping factor cf = 30 
for their models of dwarf galaxies at high redshift; Luo et 
al. (2011) adopt Cf = 20 to model high-redshift damped 
Lyman-alpha systems. 

In this paper, we adopt a variable clumping factor that 
depends on gas-phase metallicity. We use the following equa- 
tion to describe the relation between gas-phase metallicity 
and clumping factor 

Cf = Z'- (hr (15) 

for 0.01 < Z' < 1 and c f = 1 for Z' > 1. The clumping 
factor for log 10 Z 1 > — 1 is between 1 and 5, which agrees 
with the values suggested for normal galaxies in KMT09. 

In the four bottom panels of Fig. [4] we plot central 
galaxy gas-phase metallicity as a function of halo mass for 



models that adopt a variable Cf. Interestingly, the differ- 
ence between MS-I and MS-II results is largely eliminated 
for galaxies in haloes with larger than 10 11 Mq, because the 
higher clumping factor compensates the low molecular gas 
surface densities in dwarf galaxies and high redshift galax- 
ies with low gas-phase metallicities. Of course, this does not 
mean that our clumping factor prescription is correct! It just 
means that we are able to match together the results from 
the two simulations with greater ease. 



2.6 The mixing of metals 

In the models of DLB07, Guoll and FulO, metals produced 
by star formation in the disk are mixed instantaneously with 
the cold gas in the disk. SN reheating results in part of the 
cold gas (along with the metals that have been mixed into 
it) being transferred to the hot gas in the halo. If the SN 
explosion energy £sn is larger than the energy required to 
reheat the cold gas .Ereheated , the remaining energy _Esn — 
^reheated will eject part of the halo hot gas out of the halo: 
this corresponds to the so-called "ejected component" . 

We note, however, that there is direct observational ev- 
idence from deep X-ray spectral imaging of nearby star- 
forming disk galaxies that the hot gas produced by SNII 
is highly metal-enriched (Martin, Kobulnicky & Heckman 
2002) and also enhanced in alpha elements (Strickland et al. 
2004). Even in "normal" spirals, the gas is directly observed 
at scale heights 4-8 kpc above the galactic disk, before its 
emissivity decreases below the detection threshold (Strick- 
land et al. 2004) . The scenario favoured by current observa- 
tions is that SN feedback in the disks of star-forming galaxies 
create exponential atmospheres of hot gas via blow-out and 
venting of hot gas from the disk. 

We have experimented with mixing a fraction / z ,hot of 
the metals produced by star formation directly with the hot 
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Figure 5. The radial gas-phase metallicity gradient for Milky 
way sized galaxies from the model results at z = 0. The 3 dif- 
ferent colour curves represent the model results with / z hot = 
0.0,0.8, 1.0, respectively. 

gas in the halo. This turns out to have very strong impact on 
the predicted gas-phase metallicity gradients in the galaxy. 
In Fig. [5l we plot the mean radial gas-phase metallicity gra- 
dients for Milky Way- type galaxies at z = 0. The blue curve 
shows results for / z ,hot = 0.0, which is effectively the DLB07 
and Guoll mixing prescription, in which none of the metals 
are directly mixed into the hot gas. This yields a steep gradi- 
ent, with more than 0.6 dex change in metallicity in the inner 
15 kpc of the galaxy. The green curve shows what happens if 
we adopt / Zj hot = 1.0, i.e. all metals are mixed with the hot 
gas. In this case, the metallicity gradient nearly disappears, 
because we assume the hot gas in the halo is fully mixed 
and the infalling gas onto the disk has a uniform metallicity 
at different radii. According to the yield data compiled by 
Marigo et al. (2001) for a Kroupa IMF, about 20 percent 
of the total metal yield comes from AGB stars of low and 
intermediate masses. The red curve in in Fig. [5] shows the 
case where metals from AGB star outflows are retained in 
the cold ISM and metals from core-collapse SN explosions 
are mixed into hot gas (/ z ,hot = 1 — 0.2 = 0.8). 

In Sec. [5] we will compare the model gas-phase metal- 
licity gradients to observations and also explore how the pre- 
dicted metallicity profiles depend on galaxy properties such 
as mass and gas mass fraction. We will show that / z ,hot = 0.8 
gives reasonably good agreement with observations. 

2.7 The parameters in the models 

As discussed in FulO, we tune the model parameters to fit 
the observed stellar, HI and H2 mass functions at z — 0. In 
the top part of Tab. [T] we list the model parameters which 
change in this paper with respect to those listed in previ- 
ous papers. The value of H2 star formation efficiency «h 2 , 
the amplitude of SN reheating efficiency e, and the quies- 
cent hot gas black hole accretion rate k are tuned to fit the 
amplitudes of both the stellar and the gas mass functions. 
The normalization of the radial gas inflow rate a, in Sec. 
12.3.31 is tuned to reproduce the gas surface density profiles 
of Milky Way-type galaxies and to fit the observed HI and 
H2 mass functions. The fraction of metals directly mixed 
with halo hot gas is chosen to be 0.8. Other parameters are 
the same as in previous papers (bottom section of Tab.[T]): 
/b, /bh, ^merger, Y are the same as in Table 1 in Croton et 
al. (2006) (they remain unchanged in DLB07, FulO, Guoll 
and Guol3); R is same as in DLB07 (it remains unchanged 



in Guoll & Guol3); /3i, Kcheat, f], P2, Viject are same as in 
Table 2 in Guol3; 7 is the same as in Table 1 in Guoll (it 
remains unchanged in Guol3); £, Po, 01 p are the same as in 
Table 1 of FulO. 



3 STELLAR AND GAS MASS FUNCTIONS AT 
Z = 

The FulO models use the SN feedback prescriptions in 
DLB07, resulting in stellar, HI and H2 mass functions that 
are too steep at low masses (see Figure 1 in Fu et al. 2012). 
In this paper, we use the prescriptions in Guoll, which pro- 
vide a better fit to the low-mass end of the stellar mass 
function. 

In Fig.[6l we show the mass functions of stars, molecular 
gas and atomic gas for our revised models. Red solid curves 
show results for H2 fraction prescription 1 and green solid 
curves for H2 fraction prescription 2. We have computed the 
mass functions using the combined results from running the 
code on the MS and MS-II simulations. Blue circles show 
the stellar mass function derived from the Data Release 7 
of SDSS by Li & White (2009); black diamonds show SDSS 
Data Release 4 results from Baldry et al. (2008). The ob- 
served H2 mass function was derived from the FCRAO sur- 
vey (Young et al. 1995) assuming a constant CO-H2 con- 
version factor by Keres, Yun & Young (2003). The HI mass 
functions are from Zwaan et al. (2005) (blue circles) and 
Martin et al. (2010) (black diamonds). 

As mentioned in Sec. 12.71 we tune the model parame- 
ters to fit the HI, H2 and stellar mass functions, so it is no 
surprise that the agreement with models is very good. 



4 STAR FORMATION GRADIENTS 

In this section, we compare the star formation rate gradients 
in our model galaxies with observations. This provides an 
important test of whether the present-day growth of disks 
predicted by the models agrees with the data. 

We have adopted two observational samples for the 
analysis. Leroy et al. (2008) analyzed star formation gra- 
dients in a sample of 23 nearby spiral galaxies from the 
THINGS/HERACLES survey, finding that the scale length 
of the star formation rate surface density profile ?sfr is 
roughly proportional to the scale length of the stellar disk I* 
Osfr = (1±0.2)Z*). We also make use of star formation rate 
surface density measurements from the MMT long-slit ob- 
servations of 174 galaxies from Moranl2. 119 of these galax- 
ies (i.e. around 70% of the total sample) have emission lines 
strong enough to measure both SFR and gas-phase metal- 
licity gradients out to 2rso. 

In Fig. we plot the star formation rate surface den- 
sity as a function of scaled radius for model galaxies in 
four different stellar mass bins: 9.0 < log lo [M*/M ] < 9.5, 
9.5 < log lo [M*/M ] < 10.0, 10.0 < log 10 [M*/M ] < 10.5, 
and 10.5 < log lo [M*/M ] < 11.0. For the model sample, r d 
is calculated by fitting exponentials to the stellar profiles of 
the disk. The red solid and green dashed curves show results 
of mean radial profiles for H2 fraction prescriptions 1 and 
2, respectively, and the error bars on the red solid curves 
represent the ±lcr scatter between different model galaxies. 
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Table 1. The model parameters introduced/changed in this paper (top section) and other parameters whose values are the same as in 
previous papers (bottom section). Note that the parameters marked as "Croton et al. (2006)" remain unchanged in DLB07, FulO, Guoll 
and Guol3. 



Parameter 


Value 




Description 


Remark 


OCv 


0.7 km s" 


-1 T — 1 

1 kpc 1 


the ratio of radial gas inflow and gas radius 


Eq. lflO|l 


«H 2 


5.3 X 10~ 


10 —1 
yr 


H2 star formation efficiency 


Eq. Ulll 


e 


5.0 




Amplitude of SN reheating efficiency 


lab. 1 m Guoll & lab. 2 in (juolo 


"AGN 


1.5 x 10~ 


°M yr 


Quiescent hot gas black hole accretion rate 




/z,hot 


0.8 




Fraction of metal elements from quiescent star 
formation directly mixed with halo hot gas 


___ 
Sec. 12.61 


fb 


0.17 




Cosmic baryon fraction 


Croton ct al. (2006) 


/bh 


0.03 




Merger cold gas BH accretion fraction 


Croton et al. (2006) 


J merger 


0.3 




Major merger mass ratio threshold 


Croton et al. (2006) 


Y 


0.03 




Yield of metals produced per unit star formation 


Croton et al. (2006) 


R 


0.43 




Instantaneous recycled fraction of star 
formation to the cold gas 


DLB07, Guoll & Guol3 


Pi 


3.2 




Slope of SN reheating efficiency 


Guol3 


^reheat 


80 km s~ 


i 


Normalization of SN reheating efficiency 


Guol3 


V 


0.18 




Amplitude of SN ejection efficiency 


Guol3 


fa 


3.2 




Slope of SN ejection efficiency 


Guol3 


^eject 


90 km b~ 


i 


Normalization of SN ejection efficiency 


Guol3 


7 


0.3 




Ejecta reincorporation efficiency 


Guoll & Guol3 


f 


1.3 




Warm-phase correction factor 


FulO 


Po, atp 


5.93 x 10 


~ 13 Pa, 0.92 


Constant and index of the relation 


FulO 



between molecular ratio and ISM pressure 




log 10 [M,/M @ h- 2 ] log 10 [M H jM @ h- 2 } log 10 [M m /M e h- 2 } 

Figure 6. The stellar, H2 and HI mass functions from model galaxies at 2 = compared with the observations. The red solid curves 
and green solid curves in each panel are from the models with H2 fraction prescription 1 and 2, respectively. The observed stellar mass 
functions are from Li & White (2009) and from Baldry ct al. (2008). The observed H2 mass function is from the FCRAO CO survey 
by Keres ct al. (2003). The observed HI mass function is from Zwaan et al. (2005) using HIPASS data, and Martin et al. (2010) using 
ALFALFA data. 



In each panel, black curves with error bars are star for- 
mation surface density profiles from Leroy et al. (2008); Td 
for each galaxy is taken from Table 4 of that paper. The light 
blue dots are from Moranl2, who obtained long-slit spectra 
of 174 star-forming galaxies with stellar masses greater than 
10 10 A/ o from the GALEX Arecibo Sloan Digital Sky Sur- 
vey (GASS) survey (Catinella et al. 2010). In this case, 
is estimated from rgo assuming an exponential profile (i.e. 
?~90 = 3.9r<j). Note that we have simply plotted the star for- 
mation rate surface density measured for each spectral bin 
along the slit, i.e. these data points are not averaged in ra- 
dial bins as for the Leroy et al data, so the scatter from 
one galaxy to another will be substantially larger. The grey 
shaded region shows the la scatter around the median for 



the Moranl2 sample. As can be seen, the agreement between 
the Moran et al. and the Leroy et al. data is quite good. 

Our results for H2 fraction prescription 1 agree well 
with the data, particularly for the 3 highest stellar mass 
bins. H2 fraction prescription 2 yields higher star formation 
rate surface densities in the outer disks . This is because 
the pressure-based prescription yields higher H2 fractions in 
regions of low gas surface density with low gas-phase metal- 
licities. 

We note that we tuned the radial inflow prescriptions 
to match the H2 surface density profiles of the galaxies in 
the Leroy et al. sample for stellar masses log lo [Af»/M0] ~ 
10.6. Our star formation prescription is also motivated 
by results obtained for this sample. This means that star 
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Figure 7. The radial profiles of star formation surface density for model galaxies at z = compared with observations from Lcroy 
et al. (2008) (black curves with error bars) and Moranl2 (blue dots). The gray areas represent the ±lfr deviations around the median 
values for the Moranf 2 data. The red solid curves show model results of mean radial profiles for H2 prescription 1 and the green dashed 
curves are for H2 fraction prescription 2. The error bars on the red solid curves represent ±1<t scatter about the mean values for the 
model. The panels show results in 4 stellar mass bins: 10 9 - < M«/M < 10 95 , 10 95 < M*/M < 10 100 , 10 100 < M,/M Q < 10 105 , 
10 10,5 < M*/M < 10 110 . The radii plotted on the x-axis are scaled by dividing by the disk scale length r^. 



formation rate surface density profiles for galaxies with 
log 10 [M,/Mq] ~ 10.6 will match the data, essentially by 
construction. The main check in this section is to test 
whether we can match to SFR profiles over a large range 
in stellar mass. One significant conclusion that we reach is 
that H2 fraction prescription 1 does a better job than H2 
prescription 2 at matching the outer gas profiles of galaxies 
with stellar masses less than 1O 1O M . 



5 GAS-PHASE METALLICITY GRADIENTS 
IN GALAXIES 

In this section, we will compare the gas-phase metallicity 
gradients predicted by the model with data from Moranl2. 
We will first study how the gradient depends on the stel- 
lar mass of the galaxy, and then examine its dependence 
on the global gas fraction. When we compare the model 
with data from Moranl2, we choose galaxies that have HI 
and H2 gas mass fractions greater than ~ 1 — 3 %, such 
that HI and CO lines would be detected in the GASS 
and COLD GASS surveys: i.e. log 10 [M H i/M»] > -1.82 for 
galaxies with log lo [M,/M ] > 10.3, and log 10 [Af H i/M»] > 
-1.066 log 1O [M,/M ] + 9.16 for galaxies with 10.0 < 
log 1() [M,/M ] < 10.3 (Kauffmann et al. 2012). This is a 
slightly more stringent cut than adopted by Moranl2 for 
his emission-line analysis, but we have verified that the pre- 



cise location of the cut makes negligible difference to all the 
results presented in this paper. 

We adopt the gas phase oxygen abundance based on 
the 03N2 empirical index described by Pettini & Pagel 
(2004), which is derived from the N [II] A6583/He* and 
O [III] A5007/H/3 emission line ratios (see Section 3.1 of 
Moranl2 for more details). 

In the models, we only track total metallicity; the value 
of Z^e/Zq represents the cold gas- phase metallicity in units 
of the solar value. To make a comparison to oxygen abun- 
dance from observations, we adopt 12 + log 10 (O/H) = 8.69 
as the solar value (Asplund et al. 2009). We convert Z gas to 
(0/H) gas using the equation 

12 + log 10 (O/H) gas = log lo [Z gas /Z ] + 8.69 (16) 

5.1 Gas-phase metallicity gradients as a function 
of stellar mass 

In Fig. [5] we plot mean gas-phase metallicity profiles for 
model galaxies in 3 stellar mass bins: 10 9 ' 5 < M«/M < 
10 10 , 10 10 < M*/M < lO 10 ' 5 , and 10 10 ' 5 < M*/M < 
10 11 ' . The top two panels show results for H2 fraction pre- 
scription 1 and the bottom two panels are for H2 fraction 
prescription 2. In the left panels, metallicity is plotted as a 
function of radius in kpc, and in the right panels it is plot- 
ted as a function of radius scaled by rgo, the radius enclosing 
90% of the total stellar mass of the galaxy. The main result 
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Figure 8. The mean radial profiles of gas-phase metallicity from model galaxies at z = in 3 mass bins: 10 9 ' 5 < M«/M < 10 10 (blue), 
10 10 < M*/Mq < 10 10 - 5 (green), 10 10 5 < M*/M < 10 110 (red). The top two panels are the results from H 2 fraction prescription 
1, and bottom two panels are from H2 fraction prescription 2. The radii plotted on the x-axis are scaled by divising by the disk scale 
length r d . 



is that the gas radial metallicity profiles are flatter in higher 
mass galaxies; this is seen both in the left and in the right 
panels, where the radius has been scaled by the size of the 
galaxy. 

We have investigated the cause of varying metallicity 
gradients in the models by examining how the metallicity 
gradient in a central galaxy of a given halo evolves over time. 
We have found that the evolution of the gas-phase metallic- 
ity gradient is most closely tied to the merger history of the 
galaxy. The starbursts induced by the mergers consume all 
the cold gas (in the case of major mergers) or a large fraction 
of it (in the case of minor mergers) . The gas that later cools 
to form the disk has already been significantly enriched, and 
as a result, the gas-phase radial metallicity gradient in the 
new disk will be much weaker. The models thus predict that 
gas-phase metallicity gradients correlate more strongly with 
the bulge mass fraction of the galaxy than with its stellar 
mass. The correlation with stellar mass arises because mas- 
sive galaxies have experienced more mergers and have larger 
bulge mass fractions than less massive galaxies. 

This is illustrated in detail in Fig. [5] which shows re- 
sults from models with H2 fraction prescription 1. In the top 
left panel, the black solid curve shows the mean gas-phase 
metallicity gradient (in units of dex/rgo) as a function of 
galaxy stellar mass. The black dashed curves indicate the 
la scatter around the mean. As can be seen, the mean gra- 
dient increases from -0.35 dex/rgo for galaxies with stellar 
masses M, ~ 1O 9 ' 5 M to values near zero for galaxies with 
M* > 1O 1O ' 5 M . In the top middle panel, the gas-phase 



metallicity gradient is plotted as function of bulge-to-total 
(B/T) mass fraction of the galaxy. A slightly stronger, and 
more linear correlation is seen with B/T than with stellar 
mass. In the top right panel, we plot the average bulge mass 
fraction as a function of stellar mass, which shows that more 
massive galaxies tend to have higher bulge mass fractions 
(and hence weaker metallicity gradients). 

In the bottom right panel, the gas-phase metallicity gra- 
dient is plotted as a function of B/T for model galaxies in a 
narrow stellar mass interval (10.0 < log lo [Af*/M ] < 10.2) 
and in the bottom left panel metallicity gradient is plotted 
against M* in a narrow interval of B/T (0.25 < B/T < 0.3). 
These two panels demonstrate that the gas-phase metallicity 
gradient correlates primarily with bulge mass fraction. 

The purple dots superposed on the top three panels of 
Fig. [9] show data from Moranl2. We convert the concen- 
tration index R90 /R50 measured for galaxies in this sample 
to a rough estimate of B/T using the fitting equation in 
Gadotti (2009) i?go/i? 5 o = 1.93 + 2.02 B/T. Comparing the 
purple dots with black curves, the qualitative trends are in 
line with the model predictions, but the scatter is very large 
- once again, this is because the data points represent mea- 
surements along a 1-dimensional slit, rather than averages 
in radial bins as in the models. The sample is too small to 
investigate trends in metallicity gradient as a function of 
bulge fraction in a narrow stellar mass bin or as a function 
of stellar mass in a narrow bin of bulge-to-total ratio. It is 
thus not possible to assess whether the observed gas-phase 
metallicity gradient is more intrinsically correlated with stel- 
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Figure 9. Top left panel: gas-phase mctallicity (in units of dex/rgo) vs. stellar mass; top middle panel: gas-phase metallicity vs. bulge 
to total (B/T) ratio; top right panel: the B/T ratio vs stellar mass; bottom left panel: gas-phase metallicity vs. stellar mass for model 
galaxies with bulge-to-total ratio in the range 0.25 < B/T < 0.3; bottom right panel: gas-phase metallicity vs. B/T ratio for model 
galaxies in the stellar mass range 10.0 < log lo [M„/Af0] < 10.2. In each panel, the black solid curves show the mean values from the 
model sample. The black dashed curves indicate the la scatter around the mean. The purple dots show the results of individual spectral 
bin measurements from the Moranl2 data set. 



lar mass of B/T. In addition, we note that there are some 
galaxies in the observational data set with positively sloped 
radial metallicity gradients. This is not seen in the models. 
Integral field spectroscopic observations will be required in 
order to ascertain whether the outliers with positive gradi- 
ents are still present in similar numbers once metallicity is 
averaged radially. 

5.2 Relations between ^ gas , sSFR and /x„ in the 
inner and the outer regions of galaxies 

In this section, we will examine correlations between gas- 
phase metallicity, specific star formation rate and stellar 
mass surface density. In particular, we will ask whether these 
relations are the same or are different in the inner and the 
outer regions of galaxies. In a simple "closed-box" model, 
where gas is transformed into stars at a rate regulated by 
its density, chemical enrichment proceeds in the same way 
in the inner and outer regions of the galaxy. The main dif- 
ference is that less gas is consumed into stars in the outer, 
lower density regions, resulting in lower gas-phase metallici- 
ties. In our semi-analytic models, where the inner regions of 
disks assemble before the outer regions, and where gas in- 
flows and outflows regulate gas content and metallicity, star 
formation/metallicity correlations should be very different 
in the inner and outer regions of galaxies. 

Moranl2 examined relations between gas-phase metal- 
licity versus specific star formation rate and stellar sur- 



face mass density for spectral bins located in the inner 
(R < 0.7i?9o) and the outer (R > 0.7-Rgo) regions of galax- 
ies. It is difficult to draw clear conclusions from their anal- 
ysis, because the individual spectral bin measurements in 
the outer regions exhibit so much scatter. In this section, 
we work with SFR-weighted gas-phase metallicities aver- 
aged over the inner and the outer disk, which yield con- 
siderably less scattered results. We note that the individual 
gas-phase metallicity measurements are well-constrained, so 
the scatter is a real physical effect, likely to do with longer 
timescales for mixing of metals in the outer disk. In our 
models, heavy elements injected into one radial bin are in- 
stantaneously mixed throughout that radial bin, so working 
with quantities that are averaged over a larger region of the 
disk should yield a fairer comparison. We restrict our anal- 
ysis to models that use H2 fraction prescription 1 in this 
section (very similar results are obtained for the other pre- 
scription). The inner disk region is defined as r < 0.7rgo, 
and the outer disk region as r > 0.7r9o. 

In the top panels of Fig. 1101 we plot the relation between 
SFR-weighted gas-phase metallicity and mean specific star 
formation rate (sSFR=SFR/M„) in the inner disk (left) and 
in the outer disk (right). The blue dots show results for 
the Moranl2 sample and the red and yellow contours show 
results for the models. We have taken care to compute the 
average inner and outer metallicities and sSFR in the same 
way in the models as in the data. 

It is now clear that in the Moranl2 sample, the 
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Figure 10. The relation of SFR-weighted mean metallicity vs. mean stellar mass-weighted stellar surface density fi t and mean specific 
star formation rate SFR/M* for inner and outer disks. The left two panels are the results for inner disks and right two panels are for 
outer disks. The blue dots are from Moranl2 data and the contours indicate the fraction of model galaxies located in a given region of 
parameter space, as given by the colour key at the top of the plot. 



metallicity-sSFR relations are quite different in the inner 
and the outer regions of the galaxy. At fixed sSFR, the gas- 
phase metallicity is systematically higher in the inner disk 
compared to the outer disk. The systematic offset increases 
towards higher values of SFR/M,. In the inner region of 
the galaxy, gas-phase metallicity and SFR/M* are weakly 
correlated. Regions with the highest specific star formation 
rates have slightly higher gas-phase metallicities. However, 
the opposite is true in the outer regions of the galaxy: gas- 
phase metallicity decreases at higher specific star formation 
rates. The data and the models agree quite well. Q 

The reason for the "inverted" metallicity-sSFR relation 
in the outer disk can be clarified by examining the bottom 
two panels of Fig. 1101 where we plot gas-phase metallicity 
as a function of mean stellar mass-weighted stellar surface 
mass density in the inner (left) and outer (right) regions of 
galaxies for both models and data. Because stellar surface 
density is a very steeply declining function of radius in disks 
(top left panels in Figure 2), the range of stellar surface den- 
sities in inner and outer disks is almost completely disjoint. 
In the data, gas-phase metallicity correlates with stellar sur- 

1 We note that the sSFR can no longer be estimated at all accu- 
rately below a value of ~ —12, so the extension of the data points 
to sSFR values of ~ —13 should not be regarded as a significant 
discrepancy. 



face mass density only in the outer regions of galaxies, where 
stellar surface densities are low. In the models, there is a cor- 
relation between metallicity and /x, in both the inner and 
outer regions, but it is nevertheless very much stronger in 
the outer disk. Our "inside-out" disk formation models pre- 
dict that most of the stellar mass in the inner disk formed at 
high redshifts by cooling and collapse of gas within a denser 
progenitor halo. The outer disks are still in the process of 
formation at present as the halo continues to accrete dark 
matter and higher angular momentum gas is able to cool. 
The outer disks with the lowest stellar surface mass densi- 
ties are found in galaxies where gas has been accreted, but 
has not yet reached high enough densities to form molecules 
and stars. The fact that the outer gas-phase metallicities in 
models do not reach metallicities much below 0.4 solar arises 
because gas that cools from the surrounding halo has been 
significantly pre-enriched with heavy elements (see discus- 
sion in section 6). 

Finally, we would like to comment on the result in 
Moranl2 that that around 10% of disk galaxies exhibit outer 
"metallicity drops", i.e. they exhibit a sharp turndown in 
metallicity beyond ~ r$o. The magnitude of this drop was 
found to be strongly correlated with total Mht/M* ratio, but 
not with total Mh 2 /M* . We have already commented on the 
large scatter in metallicity from one spectral bin to another 
in the outer regions of galaxies in the Moranl2 sample. In- 
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deed, close examination of individual profiles in Figure 10 of 
Moranl2 reveals that there is significant diversity in profile 
shape, even among the "metal-drop" objects. 

In the top two panels Fig. 1111 we examine how the in- 
ner/outer disk metallicity difference [0/H]i n — [0/H] ou t cor- 
relates with total Mhi/M* and Mh 2 /M, ratios. Interest- 
ingly, for both the Moranl2 data and for our model galax- 
ies the metallicity difference does not correlate with either 
Mhi/M* or Mh 2 /M*. This suggests that outer metallic- 
ity drops reflect the presence of isolated spectral bins with 
low metallicity. Once again integral field spectroscopy data 
would be extremely useful to test this conjecture in more de- 
tail. Inspired by the relation between gas-phase metallicity 
radial gradient vs B/T fraction in Sec. 15. 11 we plot metallic- 
ity difference as a function of B/T in the bottom left panel 
of Fig. 1111 This anti-correlation is by far the strongest one 
in both models and data. 



6 DISCUSSION 

In this paper, we describe how we have transplanted the 
FulO prescriptions for modelling the gas and stellar profiles 
of disk galaxies and for tracking the conversion of atomic 
to molecular gas to the semi-analytic model framework of 
Guoll. Our models are run on dark matter halo merging 
trees constructed from both the MS and MS-II. Each galaxy 
disk is divided into a series of radial concentric rings. This 



allows us to track the radial distribution of gas, stars and 
metals in each galaxy. 

The main changes with respect to the FulO recipes are: 

(i) We adopt a simple star formation law in which the star 
formation rate surface density is proportional to the molec- 
ular gas surface density Esfr oc £h 2 , rather than the two 
regime star formation model in FulO. 

(ii) In Krumholz et al. H2 prescription, we adopt a 
metallicity-dependent gas clumping factor so that the gas 
is assumed to be more clumpy in low metallicity galaxy re- 
gions. This accelerates molecule formation, star formation 
and metal production in low mass galaxies and in low metal- 
licity regions and allows us to obtain convergent results in 
simulations with different resolutions. 

(iii) SNe feedback processes are no longer assumed to be less 
efficient in the regions with higher gas surface density. In- 
stead, we include radial gas inflow so that cold gas from the 
outer disk moves inwards to compensate the gas consump- 
tion in the inner disk. 

(iv) The prescription for the mixing of heavy elements pro- 
duced by star formation has been modified. Instead of mix- 
ing all the metals directly with the cold gas in the disk, we 
mix 80 % of the metals with the hot gas in the halo, and 20% 
of the metals with the cold gas in the disk. This partition 
corresponds roughly to the fraction of total metals produced 
by SNe as compared to AGB stars. 

It is worth emphasizing which of these changes were 
required in order to fit observational data: 
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Figure 12. The mean hot gas metallicity relative to the solar 
value as a function of the stellar mass of the central galaxy. 

Short gas consumption times in the inner disks of galax- 
ies were already problem in the FulO models. In that paper, 
we invoked gas surface density dependent SN feedback ef- 
ficiencies to prevent gas from being consumed too quickly 
in the inner disk. This is no longer a viable solution, be- 
cause when we apply such a model to the MS II simulation, 
we are unable to obtain a good fit to the low mass end of 
the galaxy stellar mass function. This agrees with the con- 
clusions of Guoll, who found that SN feedback efficiencies 
were required to be very high in low mass haloes to fit the 
galaxy stellar mass function at stellar masses less than few 
xlO 9 M . 

We note that the exact treatment we adopt for the ra- 
dial inflow of gas has only weak influence on the gas and 
stellar metallicity gradients in the outer regions of galax- 
ies. Outer gas-phase metallicity gradients are affected much 
more strongly by the fraction of metals from dying stars that 
is directly injected into the halo gas, rather than the inter- 
stellar medium of the galaxy. Metals ejected from a galaxy 
at early epochs are later re-accreted and this leads to flat- 
ter present-day gas-phase metallicity gradients. We demon- 
strate that a prescription in which 80 % of all the metals 
are injected into the halo gas provides the best fit to the 
relatively shallow observed metallicity gradients of galax- 
ies with stellar masses greater than 10 10 Mq (Kewley et al. 
2010; Rupke et al. 2010; Moranl2). We also show that such 
a prescription results in a good fit to the relation between 
gas-phase metallicity and specific star formation rate in the 
outer parts of galactic disks, which are still being built by 
gas accretion at the present day. 

In Fig. 1121 we plot the metallicity of the hot gas sur- 
rounding present-day central galaxies as a function of their 
stellar mass. As can be seen, the hot gas metallicity is pre- 
dicted to increase from around 0.1 solar for central galaxies 
with stellar masses of ~ 10 9 Mq to around 0.4 solar for galax- 
ies with masses comparable to or greater than that of the 
Milky Way. We note that significant fraction of the hot gas 
around central galaxies originates from material that has 
been ejected out of galaxies by supernova at early epochs. 
In low mass haloes, SN eject most of the gas out of the 
halo into the so-called "ejected" component (see Figure 5 of 
Weinmann et al. 2010). Re-incorporation of this ejected gas 
into the halo occurs after a few dynamical times for a Milky 
Way type galaxy, but substantially longer in dwarf systems 
(see Equation 23 in Guoll). 

The other main process that determines the strength 



of the gas-phase metallicity gradient in a galaxy is bulge 
formation through mergers. This is because most of the gas 
in the galaxy is consumed when the bulge is formed, and 
the metallicity gradient is re-established once new gas is 
able to accrete. In the models, galaxies with the strongest 
gas-phase metallicity gradients are those that have accreted 
gas in an undisturbed way over the age of the Universe and 
that have low bulge mass fractions. We have re-examined 
metallicity gradient trends in the Moranl2 sample and we do 
indeed find that the gas-phase metallicity gradient correlates 
more strongly with bulge-to-disk ratio than with any other 
property. 

We caution that we have used a very simple model 
Eq. (|10|l to describe radial gas inflows and our models also 
neglect radial transport of gas due to bar instabilities in 
disks. Some observational constraints on inflows do exist (e.g 
Levine et al. 2006; Zhang & Buta 2012), but the measure- 
ments are often complicated by non-radial motions in indi- 
vidual galaxies, and average inflow measurements for galaxy 
samples selected by stellar mass, size, presence or absence 
of bars etc. are still lacking. In principle, the stellar and 
gas-phase metallicity profiles in the inner regions of galaxies 
do provide additional constraints on the inflow prescriptions 
(see Figure 3) and the hope is that future data sets will mo- 
tivate more detailed work on this topic. 

Finally, we note that the evolution of radial gas, SFR 
and metallicity gradients to higher redshifts is a topic that 
is not addressed in this paper and that may be of interest 
to explore in future work. 
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